The characteristic analysis of phase-controlled array thermo-acoustic emission with multiple emitting surfaces

Thermo-acoustic (TA) ultrasound, particularly when combined with phased-controlled array technology, has garnered significant interest in the past decade due to its numerous advantages. This paper establishes a theoretical expression for thermo-acoustic phased array (TAPA) emission to investigate different acoustic fields based on input heat flow frequencies, quantities and distances of TA emission surfaces, area of emission surfaces, and phase changes between emission surfaces. The study finds that a TAPA with two emitting surfaces in a line pattern produces a consistent acoustic field compared to a single emitting surface arranged in a semicircle. Additionally, applying different phases on the surfaces narrows the scanning range with an increase in frequency, area of the TA emission surface, and the amount of emission surfaces, while enhancing the directivity of the TA wave. Moreover, increasing the distance between emitting surfaces in a square-shaped TAPA does not affect the ultrasound pressure of the main TA ultrasound but increases the quantity and size of side lobes. Furthermore, enlarging the area of emitting surfaces enhances the directivity of the TA ultrasound, necessitating a reduction in the distance between emitting surfaces or an increase in the area of the emitting surfaces in a square-shaped TAPA to enhance directivity. This paper provides a comprehensive study of TAPA to aid further research in this field.

spherical focusing TA emitter in detail.Hirota et al. 18 investigated the emission characteristics of TA emission from one-dimensional porous silicon arrays and showed that the directivity of TA waves can be achieved by changing the frequency of the input current.Gelloz et al. 19 arranged nine porous silicon TA emission surfaces into a 3 × 3 array to study the effect of changes in emission frequency and phase on the directivity of its TA emission, and concluded that proper adjustment of the size of the array element, the density of the array, as well as the relationship between the input frequency and the phase can be realised by concentrating and adjusting the direction of the emission density.Watabe et al. 20,21 investigated the changes in the characteristics of TA waves emitted by thermally induced ultrasound of porous silicon under pulse driving and found that the directivity of the TA pressure increases with the decrease of the pulse width, whereas the amplitude of the TA pressure is independent of the width of the pulse.Tsubaki et al. 22 found that the images obtained by thermal ultrasound detection have higher spatial resolution compared to conventional piezoelectric ultrasound.Zhou et al. 23 concludes that the CNT thin film spherical array panel can generate focused ultrasound beam with adjustable focal region, satisfactory ultrasound intensity.
However, existing research on TAPA is still incomplete, particularly regarding the effect of heat flow frequencies, number, distance, and phase changes of TA emitting surfaces.This work aims to study thermally induced ultrasonic emission from phased arrays using a thermally-mechanically coupled TA model.By deriving a comprehensive solution for the acoustic field, systematic investigations on TAPA can be conducted, offering significant insights into understanding TAPA features and guiding further experimental searches.

Theoretical formulas and solution
The derived equation of the TA emission from a surface source consists of the backing, ample and heating film, and considering the element surface i , phase difference θ i and the gas heat distribution factor, (i = 1, 2, 3…) 14,24 : where subscript b, s stand for backing and sample respectively.L s is the thickness of the sample, and r is the distance from TA source to the measure point.v g = γ P 0 ρ g is the sound velocity in gas [m/s], γ is the specific heat ratio of gas, κ i is the thermal conductivity [W m −1 K −1 ], and thermal diffusivity α i = κ i /ρ i c p,i [m 2 s −1 ], density ρ i [kg m −3 ], specific heat c p,i [J/(kg K)], thermal effusivity e i = k i ρ i c p,i , subscript i can take s, g, and f for the sample, gas, and heating film, respectively.C f = ρ f c f δ f , δ f is the thickness of the heating film, q 0 is the amplifica- tion [W/m 2 ].
In the following analysis of the characteristic of the TAPA from the single TA emission surface, the line pattern TA type and the square shape TA type, we assumed that the nanothermophone is placed in the air at normal atmospheric pressure and room temperature.Its schematic is shown in Fig. 1.

Results and discussions
The TA emission of a single emitting surface The comparison between theoretical and experiments results has been verified by Hu 24,25 .Thus, the related studies based on Eq. (1) can be directly utilized, as shown in Fig. 1.In our study, we used a three-layer structure (1) with silicon as the substrate, porous silicon as the sample, and coated with a heated aluminum film.The relevant parameters are listed in Table 1.
Using the formulas, we studied the acoustic properties of a single TA emission surface, analyzing the ratio of the lengths of different TA emission surfaces and the wavelength of the input frequency.The applied heat flow is q 0 = 25W/cm 2 , thus, the TA signals on a semicircular circumference t = 60 mm from the centre of the emitting surface is displayed.
The directivity of the TA sound is shown in the Fig. 2a respectively for a single square TA emitter with the slide length 2 mm, and the input frequency is approximately 50 kHz.In this case, the TA signal is uniformly distributed on a semicircle 60 mm from the centre of the TA emission surface, regarded as a point source.As the ratio of a/λ increases, the distribution of the TA signal in the 60 mm semicircle becomes similar to an ellipsoid under the input frequency f = 174 kHz and the a/λ = 1, which displayed in Fig. 2c.At the same time, the directivity has been revealed, and there is no interference from the side lobes.As the ratio of a/λ increased, the side lobes begin to appear, but the directivity of the TA emission of the square emitting surface is also becoming progressively evident.Additionally, the sound pressure is concentrated above the emitters, and the side lobes increased with the enlargement of the ratio a/λ.Moreover, the regularity in the above analysis also can clearly described in the Fig. 3.
The commonly-used RB model for sound directivity was compared with TA model in this section.The directivity in Fig. 2f is same to the Fig. 4a and similarly, the results in Fig. 3 is consistent with the survey in Fig. 4b.We compared these two models, and as expected, our model produces an acoustic field that matches what would be predicted by a theoretical circular piston in a rigid planar baffle, which illustrates the accuracy of the present model from the vision of acoustic field, including the directivity and side-lobs based on the single emitting surface.The sound pressure distribution of the two emitting surfaces As showed in Fig. 5a, its a simple array of two TA emission surfaces, the TA signal to the origin of the coordinates is t, and because the array is one-dimensional, the different TA sound field distribution would occur in yoz plane and the xoz plane.
We observed the distinction of the model in Fig. 5a, and the result is displayed in Fig. 6. and found that the tendency of the TA sound field distribution on the semicircular is invariable due to the increase of the emitting surface, but the amplitude of its TA sound pressure would increase, and the directivity of the TA wave is more accentuated, which was detected at t = 60 mm.
According to the Fig. 7, the tendency of TA sound pressure on the semicircular remains consistent, but the sound pressure decreases as the distance of two emitter increase from 6 to 14 mm, owing to no change in the difference in distance from the measuring point to the two emitting surfaces.www.nature.com/scientificreports/While the above investigates the semicircular sound field distribution in the yoz plane perpendicular to the line pattern formed by the emitting plane, the following investigates the semicircular sound field distribution in the xoz plane.
As depicted in Fig. 8.It was found that the distribution of the TA sound field on the semicircular became gradually complex as the heat frequency increased, and although the directivity was highlighted, the interference of the side lobes also became progressively frequent.Especially, the input frequency at f = 80 kHz, which is a reasonable frequency to enhance the directivity of the TA sound pressure.

The distribution of TA sound pressure based on the different phase
According to the model in Fig. 5b, we studied three emission surfaces, which array in a line and neglecting the distance of each other, that is d = 0, and a = 2 mm, f = 80 kHz, t = 60 mm.thus, The scanning of the TA sound field can be achieved from changing the difference of phase.

Effect of distance of emitting surfaces on TAPA
The Fig. 11 shows the behaviour of TA sound pressure with altitude angle for three TA emitting surfaces arranged in a line pattern, with a distance d between the different emitting surfaces, and on a semi-circular circumference 60 mm from the centre of the array of TA emitting surfaces, which at an input frequency f = 80 kHz.From the figure, it can be seen that increasing the distance between the TA emission surfaces would enhance the fluctuation of the TA pressure and the increase the side lobes, furthermore, the tendency for the TA pressure at the altitude angle θ = 90° to decrease, and the magnitude of the TA pressure of the side lobes on both sides increases, which decreasing the directivity of the whole phased array.By changing the distance between the emitting surfaces while keeping the phase distribution of the three TA emitting surfaces unchanged as (− 144°,0°,144°) → (− 108°,0°,108°) → (− 72°,0°,72°) → (− 36°,0°,-36°) → (0°,0°,0°) → (0°,30°,60°) → (45°,90°,135°), the same input frequency, the distance between adjacent TA emission surfaces d increased to 5 mm.The phased array of TA waves in the 0°-180° direction of the scanning is shown in Fig. 12, which illustrates that the increase in the distance of the emitting surfaces would lead to the TA sound waves of the parabolic lobes increase, and narrowing the parabolic scanning angle, and it would be interact with each other.
More lobes appear in (h) for the asymmetrical phase angle, thus, increasing the distance between the TA emission surfaces is not conducive to the scanning of TA waves in the 0°-180° direction.

The area of emission surface affects the TA waves emitting
As shown in Fig. 13, the directivity of TAPA increase clearly with the increased area, but would lead to the formation of multiple wave peaks of TA wave at the direction of 180°,and more peaks would appears as the increased area of the emitting surface.
A line pattern TAPA consisting of three TA emitting surfaces with an area of 5 mm × 5 mm is shown in Fig. 14, which shown that increasing the area of the emitting surface while keeping the frequency and spacing www.nature.com/scientificreports/ of the TA emitting surfaces unchanged, the TA wave of the phased array is more directional and the energy is more concentrated, but the range of the TA wave scanning becomes narrower and there are multiple paraflaps.

The TA waves emitting varied by different quantity of emission surfaces
The multiple TA emission surfaces are arranged closely in a line shape without spacing, with an area size of 2 mm × 2 mm, and under the condition of input frequency f = 80 kHz, by changing the number of TA emission surfaces, the TA pressure change on the half-circle circumference with the distance from the centre of the emission array surface of 60 mm is observed as shown in Fig. 15, and it can be seen from the figure that, with the increase in the number of TA emission surfaces, the entire TA phased array size also becomes larger, but the directivity of the emission array is getting better and better, then it will lead to an increase in the number of TA pressure peaks.The Figs. 16 and 17 represent a phased array of two and four square emitting surfaces with an area of 4 mm 2 arranged in a shape of a line, respectively.For a TAPA composed of two TA emitting surfaces, the phase distributions of the emitting surfaces are as follows: ).As for the TAPA composed of four TA emission surfaces, the phase distribution of each TA emission

The effect of arrangement of TA emission surfaces
We have discussed the TA emission surfaces closely arranged in a line array and researched the multiple influencing factors of this arrangement clearly.Thus, we arranged four TA emission surfaces with a square array, as shown in Fig. 18, to study the effect of the distance between the elements and the size of the area on the TA emission.In Fig. 19, while keeping the area of the individual TA emission surfaces constant, as the distance between the TA emission surfaces increased, which increasing the side lobes of the TA waves produced by this TAPA, and enlarging the distance would affect the TA pressure of the side lobes, but the main flap magnitude of the TA pressure never varies with the increase in the distance between the emission surfaces.
From the Fig. 20, it can be observed that when the square TA emission area of the unit is 1 mm 2 , there is only one wave peak on the curve of TA pressure variation with altitude angle, whereas when the area of a single square TA emission surface becomes 16 mm 2 , the curve contains three wave peaks.And increasing the area of emission   surface will grow more side lobes, but the ability of TA emission of whole phased array is more concentrated, and the directivity will be more pronounced.

Conclusion
Based on the previous theoretical study of TA emission, this chapter considers the influence of the sample and substrate.The method of point source integration is used to derive the acoustic pressure formula of the phased array TA emission, which serves as the theoretical basis for exploring the influences of different input heat flow frequencies, the quantity of TA emitting surfaces, the distance, the area, as well as the phase changes between the emitting surfaces and the arrangement of the acoustic wave distribution and scanning.The following conclusions were drawn.
(1) The TA pressure distribution on the semi-circular circumference forms a standing ellipse, no side lobes when a/λ = 1, Increasing the frequency strengthen the directivity of the TA emission but also increase the quantity of side lobes.www.nature.com/scientificreports/ (2) A phased array consisting of two emitting surfaces shows the same trend of TA distribution in the semicircle between the two emitting surfaces as a single emitting surface.Increasing the distance between the emitting surfaces decreases the magnitude of the TA pressure but does not change the trend of the acoustic pressure distribution on that semicircle.The semicircular sound pressure distribution in the vertical square-section over the centres of the two emitting surfaces becomes increasingly complex with growing frequency, and the number of side lobes increases.(3) The TA phased arrays arranged in a line pattern have the applied phase distribution in the order of (− 180°,0°, 0°,180°) → (− 144°,0°, 0°,144°) → (− 108°,0°, 0°,108°) → (− 72°,0°, 0°,72°) → (− 36°,0°, 0°,36°) → (0°,0°, 0°,0°), which becomes narrower with the increase of the frequency of the input heat flow, the increase of the area of the TA emitting surface, and the increase of the number of emitting surfaces, and would enhance the directivity of TA wave and enlarge the number of side lobes.However, increasing the distance between the TA emitting surfaces will facilitate the growing of grating lobes, which is not beneficial to the scanning of the phased array, and reducing the TA pressure at the top of the centre of the array surface.(4) For a TA phased array arranged in a square shape, the change of the distance between the TA emission surfaces never affect the sound pressure magnitude of the main flap of the TA wave.Increasing the distance only lead growing in both the number and the size of the side lobes.The increase of the area of the TA emission surface will enhance the directivity of the TA wave and the influence of the side lobes.Thus in order to enhance the directivity of a square-shaped TAPA, it is necessary to reduce the distance of the emitting surfaces or increase the area of the emitting surfaces.

Figure 1 .
Figure 1.Schematic of one TA emission surface, t is the distance from the measuring position to the midpoint of the TA emission surface line,θ is the altitude angle, a is the length of the square emission surface.

Figure 2 .
Figure 2. TA sound pressure directivity at different a/λ at the area of 2 mm × 2 mm.

Figure 3 .
Figure 3. TA sound pressure change with angle.

Figure 4 .Figure 5 .
Figure 4. TA sound pressure directivity and sound pressure based on the RB formula in the Ref 30 : (a) the directivity at 500KH; (b) sound pressure at 500 kHz when a/λ is 3.

Figure 8 .
Figure 8.Comparison of TA sound pressure change with altitude angle at different input frequencies.

Figure 11 .Figure 12 .
Figure 11.Variation of TA sound pressure with height angle θ for different TA surface distances.

Figure 13 .
Figure 13.TA sound pressure changes with the altitude angle under different area of emitting surface.

Figure 14 .
Figure 14.Scanning of TA waves in the 0°-180° direction by a phased array consisting of TA emission surfaces with an area of 25 mm 2 .

Figure 15 .
Figure 15.Variation of TA pressure with altitude angle for different numbers of TA emission surfaces.

Figure 16 .Figure 17 .
Figure 16.Scanning of a phased array of two TA emission surfaces in the 0°-180° direction.

Figure 18 .Figure 19 .Figure 20 .
Figure 18.The basic schematic of four emitting surfaces arranged in a square array.